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Neutron-skin thickness is a key parameter for a neutron-rich nucleus; however, it is difficult 
to determine. In the framework of the Lanzhou Quantum Molecular Dynamics (LQMD) model, a 
possible probe for the neutron-skin thickness (dnp) of neutron-rich “8Ca was studied in the 140A MeV 
48Ca + °Be projectile fragmentation reaction based on the parallel momentum distribution (p;/) 
of the residual fragments. A Fermi-type density distribution was employed to initiate the neutron 
density distributions in the LQMD simulations. A combined Gaussian function with different width 
parameters for the left side (Tz) and the right side (T g) in the distribution was used to describe the 
p// of the residual fragments. Taking neutron-rich sulfur isotopes as examples, ', shows a sensitive 


correlation with ônp of 48Ca, and is proposed as a probe for determining the neutron skin thickness 


of the projectile nucleus. 
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I. INTRODUCTION 


The neutron skin thickness, defined as the difference 
between the root mean square of the neutron and proton 
density distributions of a nucleus dnp =< r}? >1/? — < 
r2 >!/2 is an important parameter for neutron-rich nu- 
clei. With the advanced ability of new rare isotope facili- 
ties to produce nuclei near proton and neutron drip lines, 
a new era has commenced for nuclei with exotic struc- 
tures, particularly for neutron-rich isotopes at extremes 
B]. The exact value of ônp is important in research 
on symmetry energy [3] and neutron stars [4]. It is more 
important to study np of extreme nuclei because they 
may have a much lower density distribution in the surface 
region than stable nuclei. Along the history of study- 
ing np of neutron-rich nuclei, many probes have been 
proposed based on different theoretical frameworks and 
via different types of experiments. Projectile fragmen- 
tation reactions induced by neutron-proton asymmetric 
nuclei, as one type of heavy-ion reactions, are typical ex- 
perimental tools to investigate properties of neutron-rich 
nuclei, in which the light particle emissions [5] [6], frag- 
ment production ratios[7|[8], fragment mass or charge dis- 
tributions through information entropy analysis [9] [0], 
etc, are proposed as probes to study the neutron-skin 
of projectile nuclei to varying degrees of accuracy. Due 
to the difficulties in measuring neutron density distri- 
butions, the neutron-skin thickness of asymmetric nuclei 
are usually determined in experiments through indirect 
probes, for example, the giant dipole resonance i3, 
spin dipole resonance [I3], neutron-removal cross section 
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[L4], and parity violating electron scattering (for ?°°Pb 
[15] [16] and for 48Ca [I7]), etc. The momentum distri- 
bution of nucleons is widely used to study the structure 
and properties of atomic nuclei [[8}121]. The short-range 
correlation between nucleons can be experimentally stud- 
ied by detecting the high-momentum tail of the nucleon 
momentum distribution using bremsstrahlung y rays in 
heavy-ion nuclear reactions [22H24]. The parallel mo- 
mentum distribution (p,//) of fragments in breakup or 
few-body reactions is usually the first evidence of halo 
or skin nuclei, as in "Li [I8], PP [25], A1 [26], 3tNe 
[27]. Thus, p // of the residual fragments can also be em- 
ployed to determine ônp within the framework of optical 
models [28H30]. Experimentally, the p/, of fragments is 
usually used to determine their yields or cross-sections 
after integration [31|, which directly connects the yields 
of fragments and the nuclear density (and neutron skin 
thickness) of the projectile nucleus. This makes the p/, 
of fragment in potential be a probe for bn, of projectile 
nucleus. In this study, p;; of the residual fragments in 
projectile fragmentation reactions was investigated us- 
ing the Lanzhou quantum molecular dynamics (LQMD) 
model [32] [33]. It is suggested that the width of pj; of 
the fragment produced in the peripheral collisions is sen- 
sitive and could serve as a probe for db, of the projectile 
nucleus. 


II. MODEL DESCRIPTION 
A. The LQMD model 


The LQMD model is an isospin- and momentum- 
dependent transport model that includes all possible elas- 
tic and inelastic collision reaction channels during charge 


exchange. The temporal evolution of nucleons, hyperons, 
and mesons in a reaction system under a self-consistently 
generated mean field is governed by Hamilton’s equations 
of motion [32H35]. Based on the Skyrme interactions, 
isospin-, density-, and momentum-dependent Hamiltoni- 
ans were constructed [34]. The Hamiltonian of baryons 
consists of the relativistic energy, effective interaction po- 
tential, and momentum-related components. 


B EN \/ P? +M? + Uin + Umom, (1) 


where p; and m; denote the momentum and mass of the 
baryons, respectively. Uins comprises the Coulomb in- 
teraction and the local interaction potential. The local 
interaction potential is expressed as follows: 


Utoe = J Viele (r)]dr, (2) 


derived from the Skyrme energy density functional. 
Vioc(p) can be written as: 
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where pn, Pp and p = pn + pp are the neutron, 
proton, and total densities, respectively, and 6 = 


(Pn — Pp) / (Pn + Pp) is isospin asymmetry. The coeffi- 
cients a, B, Y, Jsur, gt52., and po were set to -215.7 MeV, 
142.4 MeV: 1.322, 23 MeV fm?, -2.7 MeV fm?, and 
0.16 fm? , respectively. E«,,(p) is the local part of the 
symmetry energy, which can be adjusted to mimic the 
predictions of the symmetry energy calculated using mi- 


croscopic or phenomenological many-body theories. 
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sym 
The values of Csym and ys are 52.5 MeV and 2.0, respec- 
tively, which correspond to hard-symmetry energy with 
a baryon density [34]. 

Umom takes the form 
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where Cr, = Cmom (1 + £), Cr, = Cmom (1 — z) (T # 
T’) and the isospin symbols T(7’) represent protons and 
neutrons, respectively. The parameters C'mom and € were 
determined by fitting the real part of the optical poten- 
tial as a function of the incident energy from the proton- 
nucleus elastic scattering data, and the obtained values 
of Cmom and e were 1.76 MeV and 500 c?/GeV?, respec- 
tively. Thus, the effective mass of the nuclear medium 
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FIG. 1. (Color online) The proton and neutron physical den- 
sity distributions in the initialization of Ca in the LQMD 
simulation with different dnp». In panels (a), (b), (c), and (d), 
Onp are 0.111, 0.120, 0.137, and 0.186 fm, respectively. 


at saturation density is m*/m=0.75. The parameter x is 
the strength of the isospin splitting, for which a value of 
-0.65 is adopted in this study, and the mass splitting is 
m;, > m% in the nuclear medium [86]. 

During the initialization of the projectile nucleus 
in LQMD, the initial coordinates of the nucleons are 
obtained by random sampling according to the two- 
parameter Fermi-type density distribution [38]: 


0 
Pi „i= n,p, (6) 


pi(r) = moe 
1+ exp (7555) 


where p? is a normalization constant that ensures that 
the integrated density distribution is equal to the num- 
ber of protons or neutrons, f; is a parameter for adjust- 
ing the diffuseness parameter [7] [89], t; is the diffuseness 
parameter, C; is the half-density radius. Nuclei with re- 
liable stability and an expected neutron skin thickness 
were selected as candidates for collisions, and the values 
of dnp were 0.111, 0.120, 0.136, 0.168, and 0.186 fm for 
the corresponding initial **Ca nuclei after LQMD initial- 
ization (see Fig. (i). The fragments were analyzed in 
phase space at t = 300 fm/c in the LQMD simulation, 
and nucleons with a relative distance smaller than 3.0 fm 
and a relative momentum smaller than 200 MeV/c were 
considered in the coalescence model [33]. 


B. Parallel momentum distribution 


The p;; of the fragments produced in projectile frag- 
mentation reactions exhibits a nonsymmetric distribu- 
tion in experiments, which can be fitted by a combined 


Gaussian function [81], 


= 2 í 
dv S x exp oor , ifp < po, (7) 
dp S x exp — eee , ifp> po. 


In the combined Gaussian function, the left and right 
halves have the same peak position, but different widths. 
S is the normalization factor, po is the peak position of 
the distribution, Tz and I R denote the widths of the left 
and right sides of the combined Gaussian distribution, 
respectively. 


III. RESULTS AND DISCUSSION 


p; of a fragment is generally influenced by the incident 
energy of the reaction, nuclear density of the reaction sys- 
tem, impact parameter, and the fragment itself (such as 
isospin). The experimental np determined from the pro- 
ton elastic scattering for 4Ca was —0.0107) 053 fm [37], 
and for *°Ca, it was 0.121+0.026(exp)+0.024(model) fm 
from the parity-violating method [I7]. To compare the 
sensitivity of p;; to the projectile nuclear density distri- 
bution, the 33P produced in the 140A MeV Ca + °Be 
(Onp = — 0.030 fm of Ca) and “8Ca + °Be (np = 0.168 
fm of *°Ca) reactions are simulated and fitted by Eq. 
(7), as shown in Fig. p;; for P in the symmetric 
40 Ca-induced reaction is also symmetric in Fz and Ip, 
whereas Tz is larger than Iz for 383P in the asymmetric 
48 Ca-induced reaction, indicating a neutron-skin effect 
for the neutron-rich projectile. 

With the varying impact parameters in projectile frag- 
mentation reactions, the p;; of a specific fragment is also 
influenced by the significant change in nuclear density 
from central to peripheral collisions [23], and the pj; of 
a fragment thus carries significant information on reac- 
tions, such as nuclear density distribution and nucleon- 
nucleon cross section [I]. The p;, values of the neutron- 
proton symmetric fragment ?4Mg produced in the 140A 
MeV “8Ca + °Be reaction within b = [0-2], [2-4], [4-6], 
and [6-8] fm are plotted in Fig. From the central 
to the peripheral collisions, p/; for 4Mg shifts from the 
low-momentum side to the high-momentum side as b in- 
creases. Both Iz and Ip depend on impact parameters. 
For different p,//, the value of Iz tends to decrease with 
an increase in b. Tz > Te (Fig. B), with the uncertainty 
of the I'r larger than that of Iz, (particularly for the 
central collisions). T'R remains constant when b < 6 fm 
whereas I’; decreases with increasing b. The reason for 
Tz > Tp may be the difference in the centrality in the 
reaction; that is, the larger the centrality, the deeper the 
collision with the target, and the significant p;; lost, as 
well as the wider width of Tz. 

Finally, the p;; values for the isotopic fragments with 
different isospins were studied to determine whether they 
are sensitive to the neutron skin thickness of the projec- 
tile. The neutron skin of the asymmetric nucleus mainly 
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FIG. 2. (Color online) The p,, of *°P fragments produced in 
the LQMD model simulated 140A MeV 4°48Ca + °Be reac- 


tions within b = [1-8] fm. The circles and triangles denote 
results for the °°P fragments in the **Ca (dnp = 0.168 fm) 
and *°Ca (énp = — 0.030 fm) induced reactions, respectively. 


The lines with different types denote the fitting results by a 
function according to Eq. (7). 


influences the products of peripheral collisions. The p// 
of the fragments was simulated for peripheral collisions 
of 140A MeV *8Ca + °Be within b = [6-8] fm. To com- 
pare p// values for isotopic fragments with different mass 
numbers, p;; per nucleon (p///A) was chosen. Based on 
the different values of ônp for “*Ca, as shown in Fig. 
the p;;/A distributions for neutron-rich sulfur isotopes 
(from °3S to 36S) are plotted in Fig. [4] together with the 
fitting results by Eq. (7). The values of Py, for p;; were 
obtained and the Iz correlations for the sulfur isotopes 
with dnp of 48Ca are plotted in Fig. |5| in which the corre- 
lations are fitted using the decaying exponential function. 
A strong exponential dependence of Iz for neutron-rich 
sulfur isotopes on 6,» of the projectile nucleus *%Ca is 
shown, which indicates that Tz can be an effective probe 
for the neutron skin thickness of the projectile nucleus 
in the projectile fragmentation reaction. Further simu- 
lations of the fragment de-excitation were performed us- 
ing the GEMINI code [40]. Compared with the obvious 
correlation between Iz and ônp in the LQMD simula- 
tions, Tz becomes constant as dn, varies in the LQMD 
+ GEMINI simulations, indicating that the GEMINI de- 
excitation erases the correlation between Ty and dnp. It 
is also noted that the correlation between Tz and ônp is 
an indirect probe for the neutron-skin thickness of the 
projectile nucleus because fragments are obtained after 
the compression-expansion process of the colliding sys- 
tem. Further investigations are also needed to study the 
Ty, ~ dnp correlation, including the de-excitation effects, 
since they may modify the intermediate mass fragments 
in projectile fragmentation reactions, as has been found 
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FIG. 3. (Color online) The p;; for *Mg produced in the 140A 
MeV “8Ca + °Be ( dnp = 0.111 fm of *8Ca) projectile frag- 
mentation reactions at different ranges of impact parameter 
from central collisions to peripheral collisions. The ranges of 
b are within [0-2], [2-4], [4-6], and [6-8] fm, respectively. The 
lines of different types denote the fitting results by a function 
according to Eq. (7). The inserted figure are for rz and [Rr 
for the p/;. 


in Refs. [7] [8]. 

Based on the eikonal distorted-wave impulse approx- 
imation (DWIA) explanation by Ogata et al. [80], the 
high-momentum side reflects the phase-volume effect ow- 
ing to energy and momentum conservation, whereas the 
low-momentum side reflects the momentum shift caused 
by the attractive potential of the residual nucleus when 
the incident energy is not very high (below 200A MeV). 
Because the impact parameters are restricted to b = [6- 
8] fm, the width of Ty reflects the Heisenberg uncer- 
tainty principle in quantum mechanism; that is, from 
Az: Ap > È, Az is inversely correlated to Ap [I]. 
Thus, it can be considered that the larger the Tz, the 
closer the nucleons contained in the fragment are to the 
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IV. SUMMARY 


In summary, a possible probe for the neutron skin 
thickness of a neutron-rich projectile nucleus was studied 
by simulating the 140A MeV 48Ca + Be reaction in the 
framework of the LQMD model. The neutron skin thick- 
ness of “°Ca was adjusted by varying the diffuseness of 
the neutron density distributions. A combined Gaussian 
function with different widths of the left (Tz) and right 
(Tr) halves was employed to fit p;; of the fragments. 
The p;; values of the fragments are influenced by the 
projectile nucleus, impact parameters, and the isospin of 
isotope. It was found that Iz of the p,, of the projectile- 
like fragments produced in peripheral collisions was sen- 
sitive to dnp, of the projectile nucleus. Considering that 
p; of fragments is easy to measure in experiments, the 
correlation between Iz of the projectile-like fragments 
and ônp of the projectile nucleus potentially provides a 
new probe for the neutron skin thickness of neutron-rich 
projectile nuclei. 
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FIG. 4. (Color online) The LQMD simulated p,//A in peripheral collisions (within b = [6-8] fm) for the 1404 MeV *8Ca + 
°Be reaction. The values of ônp for “8Ca are 0.111, 0.120, 0.136, 0.168, and 0.186 fm, respectively. Panel (a) is for °°5, (b) for 
345, (c) for S, and (d) for 3°S. The lines denote the fitting results by Eq. (7). 
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